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Abstract 

Kuznetsov et al. [d] and Kuznetsov and Pardo in [5| introduced the family of 
Hypergeometric Levy processes. They appear naturally in the study of fluctuations of 
stable processes when one analyses stable processes through the theory of positive self¬ 
similar Markov processes. Hypergeometric Levy processes are defined through their 
characteristic exponent, which, as a complex-valued function, has four independent 
parameters. Kyprianou et al. in [7] showed that the definition of a Hyper geometric 
Levy process could be taken to include a greater range of the aforesaid parameters than 
originally specified. In this short article, we push the parameter range even further. 


1 Introduction 

Recall that a (killed) general one-dimensional Levy process is a stochastic process issued 
from the origin with stationary and independent increments and almost sure right continuous 
paths. We write X = {X t : t > 0} for its trajectory and P for its law. The law P of a Levy 
process is characterized by its one-time transition probabilities. In particular there always 
exists a quadruple (g, a, a, II) where q > 0 is the killing rate, a 6 I is the linear coefficient, 
a e M is the Gaussian coefficient and II is a measure on M\{0} satisfying 

/ (1 A x 2 )II(dx) < oo, 

Jr 

which gives the rate at which jumps of different sizes arrive, such that 

E[e i2 ' Yt ] = e^ (i2) , z G M, (1) 

where the Laplace exponent i/)(z) is given by the Levy-Khintchine formula 

i>{z) = -q + az + Mv + f (e zx - 1 - zxl(\ x \ Kl )) n(dE). (2) 
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The (spatial) Wiener-Hopf factorisation of a Levy process £ with Laplace exponent -0 
consists of the equation 

00 ) = —k(—z)k(z), z G iR, 

where k and k are the Laplace exponents of subordinators H and H , respectively, this time 
in the sense that E[e^ Af/l ] = e ~ K ^ for > 0. The subordinators H and H are respectively 
known as the ascending and descending ladder height processes, and are related via a time- 
change to the running maximum and running minimum of the process £. The Wiener-Hopf 
factorisation has long been valued as the insight into a great many fluctuation identities and, 
accordingly, has proved to be enormously important in the general theory of Levy processes. 
See the books [ 6 ] and [I] for more on the factorisation and the central role it plays in the 
analysis of a rich variety of problems. 

There are relatively few instances where one finds concrete examples of Levy processes 
for which explicit identities exist for the characteristic exponent, the underlying quadruple 
(q, a, a, ff) as well as for the Wiener-Hopf factors. This is especially the case for Levy 
processes with two-sided jumps. One family of Levy processes for which this degree of 
tractability is available is the so-called hypergeometric Levy processes. As we shall see below, 
hypergeometric Levy processes are defined through their characteristic exponent, which can 
be written in terms of four gamma functions and which accommodates for four independent 
parameters. These processes were introduced by Kuznetsov et al. [I] and Kuznetsov and 
Pardo in [5] and given a slightly more inclusive definition in [7] by extending the parameter 
range. In this short article, we push the parameter range even further. 


We will now review the theory that has been developed for the hypergeometric Levy pro¬ 
cesses in |5j and |7| respectively. Note that the former reference refers to their generalisation 
of the Levy processes and in the latter reference as the extended hypergeometric class. As 
we are extending the definition of this class even further, we prefer to refer to all the Levy 
processes in this article as just hypergeometric Levy processes. 

We say that £ with parameters (0, 7 , j3, 7 ) belongs to the hypergeometric class of Levy 
processes if it has Laplace exponent which satisfies 


00 ) = 


T(1 — (3 + 7 — z) T (/3 + 7 + z) 


z G 


( 3 ) 


T(1 — P — z) T 0 + z) ’ 

In the existing literature hypergeometric class of Levy processes admits parameter combina¬ 
tions in the domain A\ U A 2 , where 

A = {/f<l, 76 ( 0 ,l)O> 0 , 7 e( 0 ,l)} 


and 


A 2 — {/3 e [1, 2], 7,7 G (0,1), 0 G [—1, 0]; 1 — 0 + /3 + 7 > 0,1— 0 + 0 + 7 > 0}. 

The authors in [Ej also showed that, in the parameter regime A \, the density of the Levy 


measure is 


r(») 


tt(x) = 


T(7-7)T(-7) 

m 

TO- 7 )T(—7)' 


, (1 / 3 +7> 2 JT 1 ( 1 +7;77;?7 _^. e x > 0, 

Sg+^)x^jr j(! + 7 , 77 ; 77 _ 7 ; e x ), x < 0, 


( 4 ) 


2 





where 2 Mi is the Hypergeometric function (which motivates the name of this class) and 

?? = l- /d + /3 + 7 + 7. 

Moreover, it was shown that the Wiener-Hopf factors of a hypergeometric Levy process are 


«(A) 


r(l-/3 + 7 + A) 
T(l-/3 + A) 


«(A) 


r(/3 + 7 + A) 
r(/3 +A) 


( 5 ) 


for 5RA > 0, and that in [7] it was shown that, in the parameter regime M 2 , whilst the Levy 
density remains the same as in (dj), the Wiener-Hopf factors are different. In that case, they 
can be written as 


«(A) 

for 5RA > 0 . 


(A-/3) 


r(l-/3 + 7 + A) 
T(2-p + A) 


«(A) 


(P ~ 1 + A) 


r(/3 + 7 + A) 
r(i + /3 +A)’ 


( 6 ) 


One can derive the form of the Levy density in (j4]) in a relatively straightforward way 
for both parameter regimes Mi and M 2 , by using the theory of meromorphic Levy processes 
(see 0 ). Recall that a Levy process belongs to the meromorphic class if its Levy measure is 
absolutely continuous with respect to Lebesgue measure and the density takes the form 

7 t{x) = l{a;>o} ^ a n p n e~ pnX + ljKO} ^2 ^nP n e priX , x eR, (7) 

n> 1 ri> 1 

where all the coefficients a n , a n , p n , p n are positive, the sequences {p n }n>i and {p n }n>i 
are strictly increasing, and p n —> +00 and p n —> +cxd as n —> + cxd . Part of the theory of 
meromorphic Levy processes relies on the fact that {£*,, ~(k'- k G N} and {p^, — Pk■ k G N} 
are, respectively, the roots and poles of the Laplace exponent ip, which necessarily interlace 
in the following sense 

• • • < -h < -C2 < -pi < -Ci < 0 < Ci < pi < C2 < P2 < ■ ■ ■ ( 8 ) 

For the case of hypergeometric Levy processes, since r(z) has simple poles at the points 
{— k: k G No}, the exponent ip in (J3]) has simple poles at the points {1 — (3 + ^/ + k, —(3 — 
7 — k: k G No} and simple roots at the points {1 — /3 + k, —f3 — k\kE No}. As an example, 
consider a hypergeometric process such that /3 = 0.9, 7 = 0.5, $ = 0.2 and 7 = 0.3. Note that 
these parameters lie in Mi. Then ip has simple poles of the form {—0.5 — k, 0.6 + k: k G No} 
and simple roots of the form {—0.2 — k, 0.1 + fc: k G N 0 }. The roots (dots) and poles (crosses) 
interlace as shown in the graph below. 
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The red roots and poles emerge from the factor 


r(i -fi + 7 -z) 

F(1 -P-z) 

and the bine roots and poles emerge from the factor 

r(ff + 7 + s) 
ns+z) 


(9) 


( 10 ) 


For the parameter choices here, the ascending and descending Wiener-Hopf factors are de¬ 
termined entirely by the red and bine roots and poles. Said another way, the two factors are 
precisely ([HD and (fIU|) respectively. 


Now consider another hypergeometric process whose parameters are /3 = 1.1, 7 = 0.5, 
j3 = — 0.2 and 7 = 0.8, i.e. belonging to the regime Ai- We now have simple poles of in fTUTi 
of the form {—0.6— k, 0A+k: k G No} and simple roots of the form {0.2— k, —0.1+fc: k G No}. 
Again, the interlacement of these roots (dots) and poles (crosses) are represented graphically 
below. Once again, the red roots and poles emerge from the factor ((9]) and the blue roots 
and poles from the factor (TTOl) . Note, however, for these parameter choices, the ascending 
and descending Wiener-Hopf factors are not determined by ([HD and (HH respectively. It is 
the positive roots and poles which determine the ascending factor and the negative roots and 
poles which determine the descending factor. I 11 this sense, the two Wiener-Hopf factors are 
given by (ED- 
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2 Extending the parameter domain further 

We first define the two sets of admissible parameters for that we would like to include in the 
definition of hypergeometric Levy processes: 

OO 

A 3 = (J{/3 G [0, 1], 7,7 G (0, l),/3 G [—(n+1), —n]|; 1-/3+/3+7+77, < 0, l-/3+^+7+7 1 > 0} 

n=0 

and 

OO 

A\ — |^J{/3 G [0, 1], 7,7 G (0, l),/3 G [ 77 , 77 . + 1]; n - /3 + /3 + 7 > 0,77 - f3 + j3 + 7 < 0} 

n=1 
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Our objective is to show that the the expression in (J3j) is a characteristic exponent of a Levy 
process when (ft, 7, ft, 7) G .A 3 U.A 4 . We shall preemptively consider two examples from these 
two parameter regimes and see how this suggests we might be able to extend the parameter 
range of hypergeometric Levy processes. 

For the first example, the parameters take values in * 4 . 3 . Fix ft = 0.8, 7 = 0.4, ft = —1.9 
and 7 = 0.1. We have simple roots of ip in (JHJ) at {1.5 — k, 0.2 + k: k G No} and simple poles 
at {1.4 — k , 0.6 + k: k G No}. Notice that we are in the case when n (in the definition of A 3 ) 
is equal to 1. The interlacement of the roots and poles still holds, see the figure below. Once 
again, we have coloured red the roots and poles that come from ([9]) and blue the roots and 
poles that come from (ITU]) . By allowing ft to be pushed further into the negative reals, we 
see root-pole pairs from (ITU]) have passed over to the positive side of the origin, but all the 
root-pole pairs from (JUJ) remain positive. Clearly (JUJ) and (TTU1) cannot represent the ascending 
and descending Wiener-Hopf factors respectively. However, in the spirit of (EJ), by trading 
linear factors between the two expressions in (J9]) and (TTU1) we can associate all of the positive 
and negative root-poles pairs to two different factors. As we shall see, this turns out to 
identify a Wiener-Hopf factorisation of the form 


r(i - ft + 7 - \z) 

(~ft - i~) 

T 

1 

M 

1 

_1 

X 

r(2 + ft + 7 + iz) 

1 

1— 1 

1 

1 

u 

T 

1 

1 

Th 

(—ft — 7 — 1 _ H)_ 

r(2 + /3 + i z) 


for zGl 







-3 -2 

1 


2 3 



For the second example, let us consider the case that our parameters take values in A4. 
Fix ft = 2.7, 7 = 0.2, ft = 0.3 and 7 = 0.5. We have now pushed ft further into the positive 
reals, which corresponds to passing root-pole pairs over to the negative side of 0 as seen 
below. In a similar way to the previous, we can trade linear factors between (JUJ) and (fTUl) to 
ensure that all the positive and negative root-pole pairs appear in different factors of ip in 
(|3]). The obvious candidate for the Wiener-Hopf factorisation would then take the form 


T(3 — ft - 1 - 7 — iz) 


F(ft + 7 + iz) (ft — 1 + iz) (ft — 2 + iz) 

F(3-ft- iz) 


r (ft + iz) (ft ~ 7 - 1 + iz) (ft ~ 7 - 2 + 1Z ) 




-5 

4 

3 -2 -1 

1 
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Theorem 2 . 1 . There exists a Levy process C, with Laplace exponent if as given in (J3]), where 
((3, 7 , (3, 7 ) take values in 

(i) A\. Its Wiener-Hopf factors are given by (|5j) and its Levy density is given by (J3J. 

(ii) A 2 ■ Its Wiener-Hopf factors are given by (El) and its Levy density is given by (j4j). 
(Hi) A 3 . Its Wiener-Hopf factors are given by 

r(l-/3 + 7 + A) TT + A) 


and 


^ T(1 - p + A) _ Af - j + A) ’ 


r(n + 1 + /3 + 7 + A) 
k(\) = - : -x- : —, A > 0. 


A > 0. 


( 11 ) 


r(n + 1 + (3 + A) 


Its Levy density is given by 


m 


0 +g)x (1 + 7 )k(v)k e 1 


7 T(X) = 


r(-7)r(77-7) 


E 

k=0 


kx 


(v-nt )k k\ 


r / _ ^ r( _^y e " (1 ” /3+7) "^i(i + 7, m v - 7; e- x ), x > 0, 

£M ~(0+i)x ( 1 + ^ )k(v)ke kx 


r(-7)r(7 - 7) 


k>n-\-1 


(V - 7 )k k\ ’ 


x < 0 . 


(iv) A4. Its Wiener-Hopf factors are given by 


«(A) = T{ A n ~ P V y t X \ A > 0 


and 


k(A) = 


r(l + n — [3 + A) 
r (/3 + 7 + A)fr (3-j + X 


n 


F(fi + A) %\P-T-i + y 


A > 0. 


Its Levy density is given by 

( T(r]) (1 -i3 + y) x (1 + 7 )k(y)k e~ kx 

r(- 7 )r( 7 - 7 ) " ( 7 - 7 )fc fc! 

r (7) ,-(i-^+ 7 )x V' (l + 7)fcWfc e~ fcx 

(V~f)k k\ 


x > 0 


7 T(x) = 


r(-7)r(7-7) 


m 


r(-7)r(r/ -7) 


E 

k>n 
n—1 

E 

/c =0 

#+i) Vl (l + ^, ); , ) _ 7 ;e») i I<0 


( 12 ) 


For all parameter sets, the process is killed at rate 

r(i - / 3 + T)r«; + 7 ) 

9 r(i - /?) r(/§) 

Furthermore, the process f has no Gaussian component. When 7 + 76 (0,1) (resp. 7 + 7 £ 
[1,2)) £ has paths of bounded variation and no drift (resp. paths of unbounded variation). 
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Remark 2.2. Note that the parameter ranges A EHG and A EHG defined in [?, Remark 3] 
are included in the parameter ranges A 3 and A 4 respectively. 


Remark 2.3. Let £ be a Levy process as defined in (12.11) whose parameters lie in A4. By 
noting that the ascending Wiener-Hopf factor belongs to the class of /3-subordinators, we 
have from standard literature (cf. [I]) and direct computation, with the help of the beta 
integral, that the Levy density of the ascending ladder height k is 


v(x) 


7 


r(i-T) 


(i 


-x^-y-l^n+l—P+rfx 


X > 0 . 


The density for its potential measure is given by 


u(x) 


r(7) e 


(n+ l-/3)x^ _ p -x^7-l 


)' 


x > 0. 


The same objects for the descending ladder exponent k are somewhat more complicated. 
Its potential cannot be written so easily in closed form, however, in principle, since 0 is a 
meromorphic Levy process, it can be written as a mixture of exponentials. 
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Proof of Theorem 


2.1 


For the proof of (i), we refer the reader to [5]. Also, the proof of (ii) can be found in [7j. We 
will now prove (iii). Our strategy will be to first show that the claimed factors are indeed 
Bernstein functions (i.e. Laplace exponents of subordinators) and that they belong together 
in a Wiener-Hopf factorisation, i.e. the function -0 in (j3J) is a characteristic exponent of a 
Levy process. 

We will first work under the assumption that l—(3+f3-\-^+n < 0 and 1—/ 3 +/ 3 + 7 +n > 0. 
Recall that the function 0 has simple poles at the points {1 — /3 + 7 + k, —/3 — 7 — k: k G No} 
and simple roots at the points {1 — (3 + k, —(3 — k: k G No}. 

For a given choice of parameters (/?, 7 ,/?, 7 ) G A 3 , there are n + 1 positive poles, say p*, 
of the form — f3 — 7 — j for j — 0,..., n with the rest of the positive poles, say p**, being of 
the form —f3 + 7 + l for l G N. Similarly, there are n + 1 positive roots, say (*, of the form 
— f3 — j for j — 0,..., n with the rest of the positive roots, say Q*, being of the form — f3 + l 
for l G N. The conditions 1 — /3 + /3 + 7<0,l — /3 + (3 + j > 0, which are included in A 3 , 
mean that the positive poles and roots interlace as follows: 


For k G N, denote the ordered (i.e. 


in order of interlacement) positive poles and roots by pk 
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and (k respectively, then, for z E C such that the left hand side is well defined, 


n 


! + Z /Ck 
1 + z/pk 


i + */cr fV 1 + */c; 

1 + z /pt f } 0 1 + z/p) 


n 

i> i 

n 1+ 

l> 1 -/3+7+i j=0 


1 + 


-|9+l 


n 


i + 


-p-j 


i + 


-p-i-j 


r(i — p + 7 + z) tt —4 ~ i + z 

r (i -P + z ) } = o -P-7-j + z 


= : k(z), 


(13) 


where ~ means “up to a constant”. Next, we need to apply Lemma 1 of [3] (see also equations 
(16) and (17) in that paper). For convenience, we reproduce it below in paraphrased form. 


Lemma 3.1. Assume that we have two increasing sequences p = {p n } n > i and, ( = {Cn}n>i 
of positive numbers, such that p n —> +oo as n —» +oo and the following interlacing condition 
is satisfied: 

Ci < Pi < (2 < P2 < ■■■ 


Define 


Then, for all z > 0, 


<t>(z) = JJ 

n> 1 


1 + z/p n 

1 + z/ C n 


Z > 0 . 


where 


fi(z) = ao(p, C) + 
1 


'y ^ a ri(Pi C)Cn 


3 Cn x 


n> 1 


fi{z) 


— 1 + zb 0 ((, p) + 


%{p, C) = lim TT —, <^(p, C) = (1 - Cn/Pn) TT 


'y ^ b n (C 1 p)pn 


e~ zx dx, 

(1 - e~ zx ) dx, 


D ~PnX 


n> 1 


n-H-oc -*■ ni, 
t=l 


k> 1 

k^n 


1 Cn/P/c 

i-wc*’ 


6«(c.rt = II vF Mc./» = -(i-^)n 


1 Pn/ Ck 


Ci ri ^+°° Cfc+i 

Moreover, a 0 (p, () > 0, b 0 ((,p) > 0 and for all n > 1 we have On(p,() > 0, b n ((,p) > 0. 


Cn ) 1 Pn/Pk 

k^n 


We see that, for A > 0, k(A) is a Bernstein function (in the notation of Lemma 13.11 
c j)(z ) = 1/k(z)) and so is A/k(A), A > 0. Recall that «(A), A > 0 is a special Bernstein 
function if and only if \/k(\), A > 0 is a Bernstein function, in which case the underlying 























Levy measure associated to either of these two is absolutely continuous with non-increasing 
density. We thus have that k(X) is a special Bernstein function. 

Now denote the negative poles by pk- For k > 1 they take the form (3 + 7 + k + n. 
Similarly, the zeros (k take the form (3 + k + n. Note that these poles and roots interlace in 
the following way 

• • • < -h < -C2 < ~pi < -Ci < 0 

Therefore 


n 


f + Z /Ck 
1 + z/p k 


n 


1 + 


_ z _ 

n+/3+7+fe 


1 + 


z 

n+/3-\-k 


^ r (n + 1 + /3 + 7 + z) 
r (n + 1 + j3 + z) 
=: k(z) 


A similar argument to the one in the previous paragraph shows similarly that k( A) is also a 
special Bernstein function. 

Next, using the relation T(x + 1) = xT(x) repeatedly, one easily verifies that, for the 
expression given in ( P ), = —k(—z)k(z). We can now appeal to Vigon’s theory of 

philanthropy p| Chapter 7] (see also Section 6.6 of [6]), which shows that ijj(z) is the Laplace 
exponent of a Levy process. Specifically Vigon’s theory of philanthropy states that if two 
subordinators have Laplace exponents </C(A), A > 0, i = 1,2, each of which has absolutely 
continuous and non-increasing Levy density, then they belong together in a Wiener-ffopf 
factorisation (they become ‘friends’). That is to say, 16 )9 £ M is the characteristic 
exponent of a Levy process. 

To show C is a meromorphic Levy process, apply [3, Theorem l(v)] in the killed case 
and p], Corollary 2] in the unkilled case. For Levy processes in the meromorphic class, it 
is known that their Levy measure has a density, n, of the form given in (JTJ) . We will now 
compute the coefficients dkPk and hkPk in the representation (J7|) and accordingly prove that 
7T is equivalent to the expression given in the statement of the theorem. 

Firstly, label the pk as follows: 

Pk = — /3 — 7 — k + 1, 1 <k < n + 1 

Pk = — /3 + 7 + fc — n— 1, k > n + 2 


i.e. the first n + 1 terms correspond to the poles that have been derived from (TTOl) and the 
second group of terms correspond to the poles that have been derived from ([9]) in the regime 
Ai. 

Using the relation 

a k pk = —Res (ij)(z) : z = p k ) 

we find that for 1 < k < n + 1 

(—l) fc_1 1 T(p — l + k) 

k ^ k ( k — 1)! T(1 — 7 — k) r (77 — 7 — 1 + k) 
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and for k > n + 2 


^kpk 


(- 1 ) 


k—n —2 


r(rj — 2 — n + k) 


(k — n — 2)! T(2 — 7 — k + n) r(?? — 2 — 7 — n + /c) 
Now consider the negative poles. For k > 1 we can write 

Pk = P + l + k + n 


and using 


a k p k = Res ('ijj(z ): z = -p k ) 

we can compute the coefficients a k p k in a similar fashion. Substituting these expressions into 
(ED an d recalling that 


2 ^ 1(1 + 7, V: V 


\ = (-*• + l)n(jj)n Z H 

{v - 7)n n! ’ 


z > 0, 


we obtain the density of the Levy measure as stated in the theorem. 


Now suppose at least one of the following two equalities holds 

l-/3 + /3 + 7 + n = 0, l-/3 + /3 + 7 + n = 0. (14) 

Then a finite number of the roots and poles are equal so cancel each other out in the prod¬ 
uct representation of the Wiener-Hopf factors. The remaining roots and poles still interlace 
in the required fashion so all of the expressions previously calculated still hold. However, 
this means that the coefficients corresponding to the poles that have vanished are no longer 
present in the calculation for n. 


One can prove part (iv) of the proposition in a similar way to part (iii) and we omit 
this proof for the sake of brevity. The expression for the killing rate follows by noting that 

q = -- 0 ( 0 ). 

The proof of the claim about bounded and unbounded variation follows along the same 
lines as the proof given in [[7, Proposition 1], however, we will say a few words about it. 
Formula 8.328.1 in |.2j implies that 

ip(iz) = 0(|z| 7+ ^), z —> 00 , z E M. (15) 

Proposition 2 from [lj then implies that £ has no Gaussian component. For the claims about 
the bounded and unbounded variation, the same proof as in El Proposition 1] still applies 
for all parameter regimes. For example, if we consider a process in A4, the finite sum given 
in the Levy density for x < 0 is of constant order, and since the density for x > 0 is a 
hypergeometric function minus this finite sum, the same asymptotic behaviour still holds for 
7 t(x). Finally, when £ has paths of bounded variation, [TJ Proposition 2] and (1131) shows that 
£ has no drift. □ 
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